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' Abstract 

> : 

^ . Results from the study of hadronic jets in hadron-hadron collisions at order 

in perturbation theory are presented. The focus is on various features of the 
internal structure of jets. The numerical results of the calculation are compared 
with data where possible and exhibit reasonable agreement. 
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Recent advances, both theoretical [|I], and experimental P], in the study of jet produc- 
tion in hadron collisions have made possible detailed comparisons of theory with experiment. 
In the Standard Model our general understanding of the high energy collisions of hadrons 
suggests that jets arise when short distance, large momentum transfer interactions gener- 
ate partons (quarks and gluons) that are widely separated in momentum space just after 
the hard collision. In a fashion that is not yet quantitatively understood in detail these 
configurations are thought to evolve into hadronic final states exhibiting collimated sprays 
of hadrons, which are called jets. These jets are then the observable signals of the short 
distance parton configurations. 

This general qualitative picture is characteristic of both perturbative QCD and the 
data. When one proceeds to a quantitative confrontation of theory and experiment, a precise 
definition of a jet must be supplied, and measured jet cross sections depend on the definition 
used. For instance, when one defines a jet as consisting of all the particles whose momenta lie 
inside a cone of radius R, then measured jet cross sections depend on R. On the theoretical 
side, in an order calculation a jet can consist of two partons instead of just one. At this 
level, then, a jet can have internal structure and jet cross sections calculated at order 
will depend on the jet definition applied. Two questions arise. First, how well does the 
dependence on the jet definition exhibited by the theoretical jet cross section match that 
of the experimental jet definition? Second, how well does the internal structure calculated 
at order compare to the internal structure of experimentally observed jets? We address 
these questions in this Letter. 

We consider the inclusive single jet cross section in pp collisions. This cross section is 
a function of the physical variables s, the total energy, E T , the transverse energy of the jet 
(E T = Esin9), and r], the pseudorapidity of the jet (77 = In cot (9/2)) and, as suggested 
above, the definition of the jet. The theoretical inclusive jet cross section also depends on 
the unphysical renormalization/factorization scale 71 and on the specific choice of parton 
distribution functions, f a /A(%, /■*)• m a Born level (a 2 s ) calculation the dependence on the 
scale jj, arises from both the parton distribution functions and the parton scattering cross 
section through the dependence of the latter on the strong coupling constant At order 

explicit factors of ln(/j) appear that serve to cancel a part of this \x dependence. If higher 
order contributions were calculated, they would tend to eliminate, successively, more of the 
71 dependence. At any fixed order in perturbation theory the residual 71 dependence acts 
as an estimator of the theoretical uncertainty associated with the truncated perturbation 
series. 

We find that the situation for the jet cross section at order o? s is a major improvement 
over the order o? s Born result both because the \i dependence is reduced and because the 
cross section exhibits a reasonable dependence on the jet definition. While the Born cross 
section exhibits monotonic dependence on /z, the higher order result is relatively insensitive 
to the value of 71 in a broad region near \i ~ Et/2. We estimate that the residual 
theoretical uncertainty, as indicated by the residual fi dependence, is ~ 10%. The uncertainty 
in the cross section due to the current uncertainty in the parton distribution functions is 
estimated]]]] to be somewhat larger, ~ 20%. In the perturbative calculations described here, 
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the effects of the long distance fragmentation processes and of the soft interactions of the 
"spectator" partons are ignored. We estimate that these uncalculated power suppressed 
effects constitute a correction of order ~ 6 GeV/E T to the cross section. Thus for jet E^s 
of order 100 GeV, as discussed here, the nonperturbative uncertainty is of the same order or 
smaller than that due to perturbative effects. 

At hadron colliders jets are typically defined^] in terms of the particles n whose mo- 
menta lie within a cone centered on the jet axis (rjj, </>j) in pseudorapidity 77 and azimuthal 

1 /9 

angle 0, [(rj n — r/j) 2 + (<f) n — 4>j) 2 ] < R. The jet angles (rjj, <pj) are the averages of the 
particles' angles, 

V.J = PT,nVn/E T ,J , 

nGcone 

0J = J2 VT,nKI 'E T ,J (1) 
n Scone 

with Et,j = J2n€concPT,n- This process is iterated so that the cone center matches the jet 
center (rjj, <pj) computed in Eq. (1). It is important to note that this jet algorithm is not yet 
fully defined since jets can overlap. In particular, it is possible for the constraints above to 
be satisfied by configurations where some particles are common to more than one cone. In 
the o? s calculation it is possible for a one parton jet to lie within the cone of a two parton jet. 
In the theoretical calculations described here|l[], we use the rule that only the two parton 
jets are included and the overlapping one parton jets are discarded. The precise definitions 
used by the various experiments differ to a greater or lesser extent from this form||. 

While the Born cross section with only a single parton per jet is R independent, the 
order o? s cross section can have 2 partons inside a jet and is R dependent. The theoretical 
expectation for the R dependence is shown in Fig. 1 along with results from CDFf||]. The 
inclusive single jet cross section is evaluated at E T = 100 GeV using the HMRS(B)|J parton 
distributions. Since dependence on R is not present in the Born cross section, this dependence 
is a lowest order result at order a^, that is, da/dR = + O{ol). One therefore expects that, 
although the cross section itself is relatively /1 independent, its slope da/dR will be quite 
strongly /i dependent. Just this behavior is indicated in Fig. 1 by the curves for the cross 
section versus R for 3 different /x values. (This figure is essentially Fig. 3 of Ref. || but with 
the correct theoretical result. The fourth, dot-dash curve and the parameter R sep will be 
explained below. The values of the corresponding .R-independent but strongly /i-dependent 
Born cross section are also indicated.) A correlated feature is that the \x dependence of the 
jet cross section changes as we vary R. While the order a\ jet cross section is relatively 
independent of \i for R « 0.7, for large R, R > 0.9, it is dominated by the order real 
emission process and becomes a monotonically decreasing function of \i much like the Born 
result. At small R, R < 0.5, the /1 dependence of the order cross section is dominated by 
the negative contribution from the virtual correction associated with the collinear singularity 
and becomes a monotonically increasing function of /i. In either regime we expect higher 
order corrections to be important so that the usefulness of fixed order perturbation theory 
is compromised. Thus we conclude that, at this order in perturbation theory, the results 
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are most stable for R pa 0.7. It is precisely this size that was used in the inclusive single jet 
cross section analysis published by CDF@. In some sense the perturbation theory is telling 
us that R = 0.7 is the "optimal size" for a jet cone, at least from the standpoint of making 
comparison with the order result. 

The comparison with the data in Fig. 1 suggests that, while the agreement between 
theory and data for R = 0.7 is quite good, the strong dependence on R exhibited by the 
data favors a small fi value, i.e., larger a s and more radiation. To make this comparison 
more quantitative we can characterize both the data and the theory curves in terms of 3 
parameters, 

a = A + B In R + CR 2 . (2) 

The parameterizations of the data and the theory for [i = Et/2 and \i = Et/4 are indicated 
in the first three rows of Table 1 (the parameter R sep will be defined below). We see that 
the theoretical value for the i?-independent A parameter is not too sensitive to fi. This is 
expected since it is a true one loop quantity, containing both the order a 2 contribution and 
contributions from real and virtual graphs at order oP s . The B and C terms, however, are 
subject to larger theoretical uncertainty since these terms express the R dependence of the 
cross section and this appears first at order a^. This is indicated by the sensitivity to /i 
found in the Table. We can naively associate the B term with correlated (approximately 
collinear) final state parton emission that is important near the jet direction and the C term 
with essentially uncorrelated initial state parton emission that is important far from the jet 
direction. 

The theoretical value of A agrees quite well with the experimental value of A. The 
agreement between the data and the fi = E T /2 theory for C is also quite good, but the 
agreement for B is worse than one would expect. This suggests that for the fi = Et/2 
theory, the amount of initial state emission far from the jet direction is about right but that 
there is not enough correlated radiation near the jet center. If we change fi to -Et/4, then 
the effective a s is larger and there is more radiation in all parts of phase space. Now B is 
larger, although still smaller than indicated by the data, while C is larger than indicated by 
the data. 

To examine this issue further and to analyze the internal structure of jets in detail, it is 
useful to consider the fractional E T profile, F(r, R, E T ) (we suppress the dependence on the 
jet direction (r)j,<f)j)). Given a sample of jets of transverse energy Et defined with a cone 
radius R, F(r, R, Et) is the average fraction of the jets' transverse energy that lies inside 
an inner cone of radius r < R (concentric with the jet defining cone). Said another way, 
the quantity 1 — F(r, R, Et) describes the fraction of Et that lies in the annulus between 
r and R. It is this latter quantity that is most easily calculated in perturbation theory as 
it avoids the collinear singularities at r = 0. Computing the Et weighted integral of the 
pp — > 3 partons +X cross section over the annulus and normalizing to -Et,j times the Born 
cross section yields the order a s contribution to 1 — F (the numerator is purely order 
while the denominator is purely order a 2 ). The result for F is plotted in Fig. 2 versus 
the inner radius r with R, — 1.0 for Et = 100 GeV and compared to preliminary CDF 
data|7J. Again curves for three choices of fi are exhibited. (The fourth, dot-dash curve 
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and the parameter R sep will be explained below.) As with the R dependence of the cross 
section discussed above, F is being calculated to lowest nontrivial order and thus exhibits 
monotonic /i dependence. While there is crude agreement between theory and experiment, 
the theory curves are systematically below the data for all interesting values of //. This 
situation suggests that the theoretical jets have too large a fraction of their E T near the edge 
of the jet (r ~ R). 

We have seen that the R dependence and the B parameter suggest the importance of 
higher orders to increase the level of associated radiation, at least near the center of the 
cone. At the same time our detailed considerations of the parameter C and of F suggest 
that the data favor a reduction of the Et fraction near the edge of the cone. Although 
these conclusions seem contradictory, there may be a consistent explanation based on a 
detailed but important physical point concerning how the jets are defined. We will present 
a preliminary discussion of this point here and present a more detailed study in a separate 
note.|| The issue is that of merging, how close in angle should two partons be in order to 
be associated as a single jet. In a real experiment such a situation is presumably realized 
as two sprays of hadrons, each with finite angular extent due to both fragmentation effects 
and real experimental angular resolution effects. If the angular separation is large enough, 
there is a valley in the Et distribution between the two sprays and experimental jet finding 
algorithms will tend to recognize this situation as two distinct jets. Recall that we expect 
for jets of Et > 100 GeV that the angular extent of fragmentation effects will be small 
compared to the defined jet cone sizes. However, the theoretical jet algorithm we are using 
will merge two partons into a single jet whenever it is mathematically possible. This includes 
the limiting configuration when two equal transverse energy partons (each with Et/2) are 
just 2R apart. The calculation counts this as a single jet of transverse energy Et with 
its cone centered between the two partons, i.e., centered on the valley. The treatment of 
this configuration in a real experiment will depend in detail on the implementation of the 
jet algorithm. To simulate the experimental algorithm in a simple way we add an extra 
constraint in our theoretical jet algorithm. When 2 partons, a and b, are separated by more 

than R sep {< 2R), R a b = [{Va - Vb) 2 + (0a _ <Pb) 2 } 1/2 > R sep , we no longer merge them into 
a single jet. The theoretical jet algorithm used above corresponds to R sep = 2R as noted 
in Table 1. As an example, the results of calculating both the R dependence and the Et 
fraction F with R sep = 1.3R and fi = Et/4 are illustrated by the dot-dash curves in Figs. 
1 and 2. Clearly the extra constraint of R sep has ensured that there is approximately the 
observed fraction of Et near the edge of the cone while the reduced \x value has increased the 
amount of associated radiation near the center of the cone and produced a larger variation 
with R. This conclusion is also verified by the last line in Table 1 where we observe that, 
compared to the first line of theoretical results, B has increased while C has remained the 
same. The values of both B and C that arise with these parameter choices are in reasonable 
agreement with the data. The jet cross section itself is relatively insensitive to the parameter 
R sep , decreasing by <10% as R sep is reduced from 2R to 1.3R with fixed /z for E T = 100 
GeV. 

In summary, the agreement between data and QCD perturbation theory at order for 
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the question of the dependence on the jet definition is a vast improvement over the situation 
that obtained at the Born level. There is good agreement between theory and experiment, 
at least for > 50 GeV and R near 0.7. On the question of the detailed structure within 
jets the qualitative agreement is good but there are important quantitative issues that seem 
to be dependent on the details of the implementation of the jet definition, especially the 
question of jet merging. Further study, both theoretical and experimental, is required to 
obtain a full understanding of this problem||. This is particularly interesting since there 
is some indication that such detailed internal jet structure can be invoked to differentiate 
quark jets from gluon jetsf|, ||. 
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was supported in part by U.S. Department of Energy grants DE-FG06-91ER-40614 and 
DE-FG06-85ER-40224. 



References 

[1] S.D. Ellis, Z. Kunszt and D.E. Soper, Phys. Rev. Lett. 64, 2121 (1990); Phys. Rev. 
D40, 2188 (1989); Phys. Rev. Lett. 62, 726 (1989); preprint OTIS 495, UW/PT-92-08, 
ETH-PT/92-22 (unpublished). 

[2] F. Aversa, M. Greco, P. Chiappetta and J. Ph. Guillet, Phys. Rev. Lett. 65, 401 (1990); 
Zeit. Phys. C46, 253 (1990); Nucl. Phys. B327, 105 (1989); Phys. Lett. 211B, 465 
(1988); Phys. Lett. 210B, 225 (1988). 

[3] CDF Collaboration, F. Abe, et al, Phys. Rev. Lett. 68, 1104 (1992). 

[4] See, for example, "Towards a Standardization of Jet Definitions", J.E. Huth et al, Pro- 
ceedings of the 1990 DPF Summer Study on High Energy Physics, Snowmass, 1990. 

[5] For a compilation of jet algorithms see B. Flaugher and K. Meier, Proceedings of the 1990 
DPF Summer Study on High Energy Physics, Snowmass, 1990. 

[6] P.N. Harriman, et al, Phys. Rev. D42, 798 (1990). 

[7] CDF Collaboration, F. Abe, et al, "A Measurement of Jet Shapes in pp Collisions at 
y/s = 1.8 TeV," preprint FERMILAB-PUB/92/167-E (1992). 

[8] S.D. Ellis, Z. Kunszt and D.E. Soper, in preparation. 

[9] J. Pumplin, "Variables for Distinguishing Between Quark Jets and Gluon Jets", Proceed- 
ings of the 1990 DPF Summer Study on High Energy Physics, Snowmass, 1990; "How to 
tell Quark Jets from Gluon Jets", preprint MSUTH91/3(1991). 



5 



Table 1: 3 parameter fits to data and calculated curves in Figs. 1,3 and 4 





A 


B 


C 




CDF data§ 




0.54 


0.28 


0.22 




= E T /2, R sep = 


- 2R 


0.52 


0.13 


0.19 


fi 


= Et/ 4, R sep = 


-- 2R 


0.47 


0.19 


0.30 


V = 


= Et/4, R sep = 


1.3R 


0.49 


0.22 


0.19 



Figure Captions 

Fig. 1: Inclusive jet cross section versus R for y/s = 1800 GeV, E T = 100 GeV and 0.1 < 
| T] | < 0.7 with jj, = Et/4:, Et/2, Et compared to data from CDF||; the dot-dash curve 
is explained in the text. Also indicated for the three /i choices are the values of the 
^-independent Born cross section. 

Fig. 2: F(r,R,E T ) versus r for R = 1.0, y/s = 1800 GeV, E T = 100 GeV and 0.1 < \r}\ < 
0.7 with /i = E T /4:, E T /2, E T compared to data from CDF|7j; the dot-dash curve is 
explained in the text. 
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